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It is known that the time evolution of a subsystem from an initial state to two later times, 
t\,t2 (t2 > ti), are both completely positive (CP) but it is shown here that in the intermediate 
times between ti and t2, in general, it need not be CP. This reveals the key to the Markov (if CP) 
and nonMarkov (if NCP) avataras of the intermediate dynamics. This is brought out based on A and 
B dynamical maps - without resorting to Master equation approach. The choice of tensor product 
form for the global initial state points towards the system-environment interaction dynamics as the 
sole cause for Markovianity/non-Markovianity. A succinct summary of the results is given in the 
form of a table. 
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There are three basic issues in open system quan- 
tum dynamics: (1) Role of initial correlations in system- 
environment state (2) Nature of the interaction between 
the system and its environment, and (3) Nature of the dy- 
namics at intermediate times. The first two issues were 
addressed in Ref. [ll-Q- The purpose of this paper is 
to show that the third issue is central to open system 
evolution. This is addressed here by making use of Choi- 
Jamiolkowski isomorphism [H', 'i'] . It should be noted that 
there is no use of master equation in this development. 

Understanding the basic nature of dynamical evolu- 
tion of a quantum system, which interacts with an in- 
accessible environment, attracts growing importance in 
recent years 0, [1] . This offers the key to achieve control 
over quantum systems - towards their applications in the 
emerging field of quantum computation and communi- 
cation While the overall system-environment state 
evolves unitarily, the dynamics governing the system is 
described by a completely positive (CP), trace preserv- 
ing map 0, @, ■ The concept of A and B dynamical 
maps was first introduced - as a quantum extension of 
classical stochastic dynamics - by Sudarshan, Mathews, 
Ran and Jordan (SMRJ) [10] nearly five decades ago. Its 
fundamental implications in open-system evolution are 
being put to experimental tests nowadays. We explore 
here the explicit structure of dynamical maps at inter- 
mediate times to bring out the Markov or non-Markov 
avataras of open system evolution. 

Markov approximation holds when the future dynam- 
ics depends only on the present state - and not on the 
history of the system i.e., memory effects are negligi- 
ble. The corresponding Markov dynamical map consti- 
tutes a trace preserving, CP, continuous one-parameter 
quantum semi-group 11, l3]. Markov dynamics govern- 
ing the evolution of the system density matrix is con- 
ventionally described by Lindblad-Gorini-Kossakowski- 
Sudarshan (LGKS) master equation [ll|, [13 ^ = Cp, 
where C is the time- independent Lindbladian operator 



generating the underlying quantum Markov semi-group. 
Generalized Markov processes are formulated in terms 
of time-dependent Lindblad generators and the asso- 
ciated trace preserving CP dynamical map is a two- 
parameter divisible map [l^, [l3| and this too corresponds 
to memory-less Markovian evolution. Last few years have 
witnessed intense efforts towards formulating, character- 



izing and quantifying non-Markovian dynamics |13l-l21 



NCP dynamical maps do make their presence felt in 
the open-system dynamics obtained from the joint uni- 
tary evolution - if the system and environment are in 
an initially quantum correlated state [H-Ql- In such 
cases, the open-system evolution turns out to be non- 
Markovian [17]. However, the source of such non- 
Markovianity could not be attributed entirely to cither 
initial system-environment correlations or their dynam- 
ical interaction. This issue gets refined if initial global 
state is in the tensor product form, in which case the 
sole cause of Markovianity/non-Markovianity is dynam- 
ics. It is this aspect that is addressed here. It may be 
noted that these issues could not be addressed in the 
generalized master equation approaches [H M, ■ 

We now begin with a brief review of the A and B dy- 
namical maps [l3] . These dynamical maps transform the 
initial system density matrix psita) to final density ma- 
trix psit) via. 



ai, 02, a[,a'2 = 1, 2, . . . , d 



where the realligned matrix B is defined by. 



B„ 



(3) 
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The identity A-map has the foUowing structure: 

^ ^aia2;a[a'2 ~ ^ai,a[ ^02, (4) 



Trace preservation : 

Hermiticity : 
Positivity : 

ai ,a2 ,02 



Clearly, the B map is positive, Hermitian (P x cP matrix 
with trace d. We would also like to point out here that the 
composition of two yl-maps, Ai * A2 is merely a matrix 
multiplication, whereas it is not so in its i?-form. 

Let us consider unitary evolution of global system- 
environment state ps(to) psito) from an initial time 
to to a final time t2 - passing through an intermediate 
instant ti (i.e., to < ti < 12)- The A-map associated with 
to to ti and that between to to t2 are readily identified 



If {U{h,to)Ps{to) ® PE{to)UHtiM)] = Psih) <E) PEih) 

(memory-less reservoir condition) then the LHS of ([7|) is 
also given by, 

Trs [C/(t2, ti)ps(ti) (g, pB(tl)[/t(t2, ti)] = A{t2,ti)ps{ti). 

(8) 

With j — 1 in ([6]), we substitute ps{to) = 
A~^{ti,to)ps{ti) in ^ to identify the intermediate A 
map as, 

A{t2,h) = A{t2,to)A-\ti,to). (9) 

We would like to emphasize that when the environment is 
passive ( Markovian dynamics) , the intermediate A- map 
has the composition as in In such cases A{ti,t2) 

is ensured to be CP - otherwise it is NCP, and hence 
non-Markovian. When pssiti) Ps{ti) ^ P£;(ti), im- 
plying system-environment correlations at any interme- 
diate time ti, the composition ^ does not hold. This 



The requirement that the density matrix ps (t) has unit 
trace, Hermitian and is positive semi-definite leads to the 
conditions 



(5) 

ai ,a2 ,0.^ ,0.2 

I 

as follows: 

Tie [U(tj,to)ps{tQ) <E} PE{to)U\tj,to)] ^ A{tj,to) ps{to) 
^psitj), J -1,2. (6) 

The dynamical map A{tj,to) is CP. 

From the composition law t/(t2, to) — ?7(t2, ti)J7(ti, to) 
of unitary evolution, we obtain. 



(7) 

I 

construction provides us with a test of the CP or NCP 
nature of the process by explicit evaluation of the RHS 
of (The corresponding intermediate i?-map B{t2,ti) 
is positive if the process is CP - otherwise it is NCP). 

We now observe that Jamiolkowski isomorphism Q 
provides an insight that the i?-map is directly related 
to a (t^ X d^ system-ancilla bipartite density matrix. 
More specifically, the action of the map A^'^ (E) A on 
the maximally entangled system-ancilla state IV'me) = 
73 'Et=o \h i) results in i B i.e., 

[A"® A] |V'me)(^me| = (10) 

gives an explicit matrix representation for the B-map 
(Here A}'^ is the identity A- map (see which leaves 
the ancilla undisturbed). 
In detail, we have. 
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ai 

/I* r> o* 

ia2:a'-^a'^ -^a^ay.a^a'-^^ ^aia^-,a2a'^ ^a2a'2-,aia'^ 

^ ^ -^ai -^0.2 ^aia2;a[a'2 Ua'^ Ua' — Oj ^ ^ •'^ai Va'^ Baj^a[;a2a'2 ^a2 UaL — ^ 
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Tl-E [U{t2,ti){U{ti,to)ps{to)<»PEito)UHti,to)}U^{t2,ti)] =A{t2, to) Ps{to)- 
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I A -Ir 

^^bib2]aia2 — ^ -'-'fci ai ;b2 ^2 • 



(11) 



In other words, Jamiolkowski isomorphism maps ev- 
ery completely positive dynamical map acting on d di- 
mensional space to a positive definite d^ x d^ bipartite 
density matrix pab (which is just ^B) - whose partial 
trace (over the first subsystem index) is a maximally dis- 
ordered state. This result is powerful - as one may now 
identify several toy models of dynamical B maps to in- 
vestigate the nature of dynamics, which will be presented 
next. 

Consider the two-qubit (system-ancilla) density ma- 



[i-p(*)] 



p{t) |'0me)(V': 



ME 



with 



trix Pab{t) = ^ 4^ -'2 12 

< p{t) < 1 being a suitable function of time, and 
IV'me) = ^ (|0a,0fc) + |la,l&)) - as a prototype of our 
dynamical B map. For a dynamical map, time depen- 
dence in p{t) occurs due to the underlying Hamiltonian 
evolution. This state is especially important in that it 
exhibits both separable and entangled states, as its char- 
acteristic parameter p{t) is varied. Its use here as a valid 
_B-map is novel in identifying Markovianity of dynamics. 

On evaluating the corresponding A map A(t,0), one 
can obtain the intermediate dynamical map ^4(^2, ^i) = 
A{t2, 0)^'^(ii, 0). The intermediate time i?-map 
B{t2,ti) is then given by 



Bit2, ti) 



[pjti) - p{t2)] . ^ . , Mh) 



|'0me)(V': 



Its eigenvalues are Ai 
A4 = p(«iH3p(*^) . A choice p{t) 



ME I- 

(12) 
and 



P(*l)-P(«2) 

\4 ~ - ■ 2p{t[) — ^ cuoice p\i) = cos^^(ai) for any 
> 1 leads to NCPness of the intermdiate map - as the 
eigenvalues Ai,2,3 of B{t2,ti) turn out to be negative 
- and hence non-Mar kovian dynamics ensues. Another 



choice p(i) 



corresponds to a CP intermediate map 



- resulting in a Markovian process. In this case, we also 
find that A{t2, ti) — A{t2 — ti) and this forms a Markov 
semigroup. However, if p(t) — e""*'', (/3 7^ 1), the inter- 
mediate map is still CP (and hence Markovian), though 
A{t2,ti) A{t2 — ti) - and hence it does not constitute 
a one-parameter semigroup. 

Another important feature that we wish to illustrate 
through this example is the following: Concurrence of 
Pab{t) = 2B{t,0) (given by C = M^ilL) can never in- 
crease as a result of Markovian evolution. This is because 
ensuing dynamics is a local CP map on the system. Any 
temporary regain of system-ancilla entanglement during 
the course of evolution is clearly attributed to the back- 
flow from environment to the system - which is a signa- 



ture of non-Markovian process. This feature is displayed 
in Fig. 1 by plotting the concurrence pab{t) for different 
choices of p{t). 
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FIG. 1: 

Pab{t) = 



1 



Concurrence C — 



2 3 4 

at 



3p(t)- 



of the system-ancilla state 
I2 <8i h +p{t) |^/'me)(i/'me1, vs scaled time at, 
for the following choices (i) Markov process: p{t) = e~"* (solid 
line) and (ii) non-Markov process: p{t) = cos^^ (at), N — 1 
(dashed line) and N = 5 (dot-dashed line). Note that there is 
a death and re-birth of entanglement (dash, dot-dashed lines) 
due of back-flow from environment. 



We now present two Hamiltonian models, which give 
rise to explicit structure of time dependence in the open 
system evolution. Interaction Hamiltonian considered 
here is 21 1 



H 



A 



N 



(13) 



k=l 



This is a simplified model of a hyperfine interaction of 
a spin- 1/2 system with N spin- 1/2 nuclear environment 
in a quantum dot. Taking the initial system-environment 
state to be ^^(O)®-^^, the dynamical A-map is obtained 

by evaluating Trs U{t,0) ps{0) ^ U^t,0) (where 
U{t,0) = Exp[-iHt]): 

A{t,0) = ^{l~x{t))a,^<j, + ^{l + x{t))l2<S)l2 



x{t) = cos^^ I 



(14) 



From this, the intermediate map A{t2jt\) (see ©) and 
in turn the corresponding B{t2,ti) may be readily eval- 



4 



Initial state 


B(0,ti) 


B{Q,t2) 


B{t2,tl) 


Nature of dynamics 


PSE = ps ® PE 


CP 


CP 


CP 


Markov 


PSE ^ ps ® Pe 


CP 


CP 


NCP 


non-Markov 


Quantum correlated psE 


NCP 


NCP 




non-Markov 



TABLE I: CP/NCP nature of the maps and the associated 
Markov/non-Markov features 



uated. We obtain, 



Bit2, h) = - {h®h+(y,®a,) 



X{t2) 

'2x{ti) 



{(JX® - (Ty ® ay) . (15) 



The eigenvalues of B{t2,ti) are 0,0, 1 ± fpYJ- Clearly, 
the intermediate time dynamics is NCP - and hence the 
process is non-Markovian. 

We now consider the open system dynamics arising 
from the unitary evolution [l[ 



U{t,0) 



i t [ui a~ ®a^\ 

cos(a;t/2) I2® I2— «sin(a; t/2) Gz ® 



(16) 



on the system-environment initial state ps£;(0) = ^5(0)® 
Ps(0) = i (^2 + Ox) ®\{J2 + Oz) ■ The A{t, 0) map is 
given by. 



^(^.0) = ^(1+C0S(^^)) /2®/2 + ^ (1 - C0S(W<)) Oz^Oz. 

(17) 

Following we obtain 



-8(^2,^1) = 2 (^2«)/2+CT2<»CT2) 
C0s(wf2) 

-K- — — {Ox ® Ox - Oy® Gy) . (18) 

2cos(a;ti) 

The eigenvalues of the i?-map being 0, 0, lifff^ bring- 
ing out the NCP nature of the intermediate map explic- 
itly - implying non-Markov nature of the dynamical pro- 
cess. This model, with initially correlated states, has 
been explored in Ref. [J [l^ and the dynamical map 
turned out be NCP throughout (not merely in the in- 
termediate time interval). 



In conclusion, a few remarks on a variety of definitions 
of non-Markovianity in the recent literature may be re- 
called here. Mainly the focus has been towards capturing 
the violation of semi-group property 16|, |17[ or more re- 
cently - its two-parameter generalization viz the divisibil- 
ity of the dynamical map Il3l . [l3 |. Yet another measure, 
where non-Markovianity [l8| is attributed to increase of 
distinguishability of any pairs of states (as a result of 
the partial back-flow of information from the environ- 
ment into the system) and is quantified in terms of trace 
distance of the states. It has been shown that the two 
different measures of non-Markovianity - one based on 
the divisibility of the dynamical map [14] and the other 



based upon the distinguishability of quantum states 18 1 
- need not agree with each other [l^ . A modified version 
of the criterion of Ref. [13] was proposed recently . In 
this paper it is shown that the intermediate time evolu- 
tion of the system is not necessarily CP and holds the 
key to the memory of the initial state in a subtle way. 
The results are summarized in Table 1. 
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